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ABSTRACT 

We present near-infrared polarimetric images of the dusty circumstellar envelope 
(CSE) of IRAS 19306+1407, acquired at the United Kingdom Infrared Telescope 
(UKIRT) using the UKIRT 1-5 |i.m Imager Spectrometer (UIST) in conjunction with 
the half-waveplate module IRP0L2. We present additional 450 and 850 |^m photome- 
try obtained with the Sub-mm Common User Bolometer Array (SCUBA) at the James 
Clerk Maxwell Telescope (JCMT), as well as archived Hubble Space Telescope (HST) 
F606W- and F814W-?i\tei images. The CSE structure in polarized flux at J- and K- 
bands shows an elongation NNE-SSW with two bright scattering shoulders NW-SE. 
These features are not perpendicular to each other and could signify a recent 'twist' 
in the outflow axis. We model the CSE using an axisymmetric light scattering (ALS) 
code to investigate the polarization produced by the CSE, and an axisymmetric radi- 
ation transport (DART) code to fit the SED. A good fit was achieved with the ALS 
and DART models using silicate grains, 0.1-0.4 |i.m with a power-law size distribution 
of a~'^'^, and an axisymmetric shell geometry with an equator-to-pole contrast of 7:1. 
The spectral type of the central star is determined to be BlI supporting previous 
suggestions that the object is an early PN. We have constrained the CSE and inter- 
stellar extinction as 2.0 and 4.2 mag respectively, and have estimated a distance of 
2.7 kpc. At this distance the stellar luminosity is ~4500 Lq and the mass of the CSE 
~0.2 Mq. We also determine that the mass loss lasted ~5300 yrs with a mass- loss 
rate of -3.4x10"'^ Mq yr'^. 

Key virords: stars: AGB and post-AGB - stars: circumstellar matter - infrared: stars 
- stars: individual {IRAS 19306-1-1407) - stars: mass loss - techniques: polarimetric. 



1 INTRODUCTION 

Post-asymptotic giant branch (post-AGB) stars are lumi- 
nous (10^-10'' Lq) evol ved stars wi th initial masses in 
the range 0.8-8 Mq (see I Van Wincke l 2003, for a general 
review). At the end of the AGB phase, mass- loss rates 
can peak at over 10"* Mq yr~^ before dropping dra- 
matically, as the star enters its post-AGB evolution (e.g. 
ISchonberneil |l983i), creating detached envelopes of gas and 
dust. These dusty circumstellar envelopes (CSEs) are then 
visible at optical and ne ar-infrared w avelengths as proto- 
planetary nebulae (PPN; |Kto3[1993). A seemingly ubiq- 
uitous feature of PPN is their lack of spherical symme- 
try, with many having a bipolar or point-symmetric struc- 
ture. Notable and well-studied exam ples are the Egg Neb- 
ula (AFGL 2688; ISahai et al.l ll^98^ and the Red Rect- 
angle (AFGL 915: ICohen et al.ll2003) . Optical and near- 

* E-mail: kloweSstar .herts . ac .uk 

t E-mail: t .gledhillSstar .herts . ac .uk 



infrared surveys of PPN have shown that in all cases where 
a CSE i s detected then it appears a symmet ric in some 
way (e.g. lUeta. Meixner fc Bobrowskvli2000: Gledhill et al.l 
12001). Possible mechanisms for the shaping of PPN usually 
involve interaction of the mass-losin g star with a b inary com- 
panion, and have been reviewed bv lBalick fc Fra nk (2002). 

Imaging polarimetry is a differential imaging technique, 
which is well-suited to the study of CSEs surrounding post- 
AGB stars. The technique discriminates between the faint 
but polarized scattered light from the PPN and any bright 
unpolarized emission from the central star. This enables the 
imaging of circumstellar material that would normally be 
lost under the wings of the stellar point spread function 
(PSF), thereby obtaining information on the dust distribu- 
tion close to the central source. Imaging polarimetric surveys 
of post-AGB stars using the UK Infrared Telescope have 
detected scattered light from PPN around 34 stars, and all 
of these PPN were found to be axi symmetric in some way 
jGledhill et al.ll200ll: lGledhil]||2005^ ■ Higher spatial resolu- 
tion polarimetry using the Hubble Space Telescope (HST) 
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has enabled more detailed studies of the morphology of PPN, 
as well as providing constraints on dust grain properties 
in these systems, and has rev ealed p oint-sy mmetries, iets 
and multi-lobed struc tures (e.g. lUet~Murakawa fc Meixnerl 
I 20051: ISu et"al]|2003ll . 

In this paper, we examine IRAS 19306-1-1407 (GLMP 
923), which has IRAS colours typic al of a post-AGB star 
with a cold CSE lOmont et al ]|1993!). Radio and millimetre 
surveys for m olecular emis sion have failed t o detect OH or 
H2O masers (iLiklffl llllQSgl) or CO emission (|Arauilla et al.l 
Il986l: iLikkel et alj Il991fl . Howeve r, the object shows a 
number of dust spectral features. iHrivnak. Volk fc Kwokl 
l|2000(l present ISO spectroscopy showing emission fea- 
tures at 6.3, 7.8 and 10.7 \i.ra, with a "probable" feature 
at 3.3 \im, and compare these features to the uniden- 
tified infrared (UIR) bands at 3.3, 6.2 and 7.7 ^.m, 
commonly attribut ed to polycyclic aromatic hydrocarbo n 
(PAH) molecules jAUamandola. Tielens fc BarkeJ Il989l) . 
Given that the mid-infrared spec tral features are simi lar to 
those seen in hot carbon-rich PN. lHrivnak et all ll200(i) sug- 
gest that the object is a young PN. A further analysis of 
the 750 data by Hodg ^ ^t a l. ( 2004) confirms the presence 
of UIR features at 3.3, 6.2, and 7.7 \im, with the addition 
of the 8.6 and 11.2 ^im features. These authors also men- 
tion the presence of silicate emission at 11, 19 and 23 \im, 
raising the possibility that IRAS 19306+1407 may have a 
mixed CSE chemistry. 

Optical spectroscopy shows a broadened Ha emission 
line with line width of ^2300 km s~^ i ndicating a fast out- 
flow llSahai fc Sanche z Contrera^|200 4r). as well as H/3 and 
[Nil] emission, leading iKelly fc HrivnaM I2OO?) to suggest 
a spectral type of approximately BO for the star. A num- 
ber of H2 emission lines are seen in the K-hanA, with line 
ratio s suggesting a mix of radiative and shock excitation 
iKel ly & Hrivnak 2005). Imaging through a narrow-band H2 
fllter, centred on the 2.122 ^m line, shows that the H2 emis- 
sion has a ring-like structure with e vidence for bipolar lobes 
exten ding perpendicular to the ring iVolk. Hrivnak fc Kwokl 
l2004l) . 

We present the first near-infrared polarimetric images 
of the dusty CSE of IRAS 19306+1407, showing the struc- 
ture of the envelope in scattered light. We also present new 
sub-millimetre photometry and archived HST images. The 
observations are interpreted using 2-dimensional (axisym- 
metric) light scattering and radiation transport models. 



2 OBSERVATIONS AND RESULTS 

2.1 Imaging polarimetry observations and results 

Polarimetric imaging at J- and K-h&nA of IRAS 19306+1407 
was obtained at the 3.8-m United Kingdom Infrared Tele- 
scope (UKIRT) on Mauna Kea, Hawai'i, using the UKIRT 
1-5 micron imager spectrometer (UIST) in conjunction with 
the infrared half-waveplate module (IRPOL2) . A pixel scale 
of 0.12 arcsec was used and observations were made on 2003 
June 8 with an average seeing of 0.5 arcsec. The total inte- 
gration time for each filter was 237.6 seconds, comprising 24 
exposures of 9.9 seconds each (see Table0. Linear polarime- 
try was obtained by observing at four half-waveplate angles 
of 0°, 22.5°, 45° and 67.5°. The data reduction was carried 



Table 1. Summary of photometry for IRAS 19306+1407 for HST 
(using Vega zero points), UKIRT and SCUBA observations, in- 
cluding integration time (Int.) and the extent (Size) of the semi- 
major and minor axes of the aperture used in photometry. The 
PA angle of photometry aperture is equal to 18° (E of N). 



Band 


Magnitude 


Flux 


Int. 


Size 






(mjy) 


(s) 


(arcsec X arcsec) 


F606W^ 


13.81 + 0.03 


9.5 + 0.3 


300 


3.2 X 2.0 


F8UW^ 


12.45 + 0.02 


26.1 + 0.5 


50 


3.2 X 1.9 




11.18 + 0.04 


51.5 + 0.8 


237.6 


3.9 X 2.4 




10.29 + 0.12 


48.4 + 2.2 


237.6 


3.9 X 2.4 


45OW 




49.9 + 38.7 


1334t 




850 W 




14.1 + 3.7 


1334t 





Notes: central wavelengths at ^0.5888|jm (Broad V), '^0.8115nm 
(Johnson 7), '=1.25|am, '^2.2\xm, =450vmi and f850nm; and 
t inclusive of observational overheads. 



out using STARLINK^ applications. A bad pixel mask was cre- 
ated using ORACDR and chopped to 512 by 512 pixels. The 
standard subtraction of dark frames and flat fielding were 
carried out by CCDPACK. A 3D cube consisting of the /, Q 
and U Stokes images, was produced using POLKA from the 
POLPACK suite, and this was then used to derive the per cent 
polarization, polarized flux and polarization angle. A more 
detailed description of dual-bea m polarimetry and the d ata 
reduction techniques is given bv lBerry fc Gledhilll Jl999D . 

Photometric standards, FS 147 (J) and FS 141 (TT), 
were used to flux calibrate the data giving J=11.18 ± 0.04 
and JC=10.29 ± 0.12. 

For these observations, the focal plane polarimetry 
mask was removed, so that a 512 by 512 pixel sub-array of 
the UIST detector could be used. This enabled faster read- 
out times and exposures of less than 1 second, so that obser- 
vations of bright sources could be made without the risk of 
saturation. This conflguration of UIST resulted in the over- 
lapping of the 0- and e-beams produced by the WoUaston 
prism and a final analysis area of 20 by 60 arcsec. The Wol- 
laston prism splits each star into an e- and o-component sep- 
arated by 20 arcsec, so that any star in the field lying more 
than 10 arcsec along the prism dispersion axis from the tar- 
get will only have one component in the analysis area. Since 
both e- and o-beams are required to correctly calculate the 
Stokes intensities /, Q and U, these offset stars appear as 
highly polarized artefacts in the reduced data, and they are 
marked as such on Fig. [H As the prism dispersion varies 
slightly with wavelength, this results in an apparent shift of 
the artefact stars between the J- and ii'-filters. 

The J- and K-hand polarimetric results are shown in 
Fig-0 The total intensity images are shown in Fig. 0(a) and 
(c), superimposed with polarization vectors, and show the 
centrally peaked nature of the source. The object is clearly 
extended, relative to the 0.5 arcsec seeing FWHM, with faint 
emission detected out to a radius of approximately 3 arcsec. 
The lowest contour in both filters is 3 times the sky noise and 
in the Jj image, shows that the faint emission is elongated in 
a NNE/SSW direction. Details of contour levels are given in 
the Figure caption. It is possible that a similar extension is 
present in the 7k image, but confusion due to the presence 
of the artefact stars makes this uncertain. 

The polarized fiux, produced by light scattering from 

^ Available from www.starllnk.ac.uk 
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Figure 1. High resolution images are available at Ihttp:/ /star-www. harts. ac.uk/~k;lowe/) The J- and if-band observations 
are displayed at the top and bottom of the figure respectively. These images have been scaled logarithmically. The total intensity (I) 
is displayed in sub-figures (a) and (c) with overlaid polarization vectors (pol). Sub-figures (a) and (c) are scaled between 20 and 13 
mag arcsec"^. The lowest outer contour levels are 19 and 18 mag arcsec"^ and separated by 1 mag arcsec"'^ for (a) and (c) respectively. 
The polarized flux (IP) images (b) and (d) are scaled between 20 to 16 mag arcsec"'^ and 19 to 16 mag arcsec"^ respectively. The lowest 
outer contours are 19 (b) and 18 (d) mag arcsec"^ and separated by 0.5 mag arcsec"^. 



dust grains, is shown in Fig. Q(b) and (d). In both filters, 
the central region appears elongated along a PA 136° East 
of North, with two bright shoulders of emission either side of 
the star. At J (IPj image) this structure is embedded within 
fainter more extended emission orientated at 18° East of 
North, seen in the lowest three contours (the lowest contour 
is at 1.5 times the sky noise). This faint extension is not as 
apparent in the K-hand polarized flux image (IPk), which 
is approximately 1 mag arcsec"^ shallower than the J-band 
data. The NW shoulder is brighter than the SE shoulder, 
particularly apparent in the IPj image. Similar morphology 
has b een observed in pola rised flux in a number of other 
PPN. iGledhiU et"ai] l|2n0ltl found bright arc-like structures 
on either side of the star in IRAS 17436+5003 as well as 
shoulder-like features in IRAS 19500-1709 and more ring-like 



features in IRAS 22223-1-4327 and 22272-1-5435. They inter- 
preted these structures in terms of scattering from the inner 
surfaces of a detached axisymmetric shell, with an equatorial 
density enhancement, and classified these objects as "shell- 
type" objects. The arcs in IRAS 17436+5003 were later fully 
resolved i n mid-infrared imaging of thermal emission from 
the dust jGledhill fc Yatesll2003ll and successfully modelled 
using an axisymmet ric dust distribution based on that of 
iKahn fc Wed (Il985ll . Further evidence for arcs and shoul- 
ders is seen in polarized flux i mages of IRA S 06530-0213, 
07430+1115 and 19374+2359 iGledhilj lioosi) and was in- 
terpreted using light-scattering in a Kahn fc West density 
distribution. We therefore interpret the polarized flux shoul- 
ders seen around IRAS 19306+1407 in the same way, and 
suggest that they result from increased scattering at the in- 
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Table 2. Summary of polarimetric results of IRAS 19306+1407 
for each band, detailing the maximum polarization, integrated 
polarization and the position angle (E of N) of the major and 
minor axis of the nebula in polarized flux. 



Band 


Max. Pol. 


Integrated Pol.T 


PA 

^ -'^major 


PA 




(per cent) 


(per cent) 


n 


(°) 


J 


15 ±6 


1.7±0.1 


18 


136 


K 


10 ±4 


1.3 ±0.1 


18 


136 



' - The integrated polarization over the source with apertures of 
radii of 1.7- and 1.4-arcsec for J- and if-band respectively. 



ner boundary of a detached shell with an equatorial dust 
density enhancement. 

The polarization vectors shown in Fig. 1 (a) and (c) are 
binned over 0.36 x 0.36 arcsec (3x3 pixels) and have a 
signal-to-noise threshold of 2 in per cent polarization. The 
vector pattern appears approximately centro-symmetric in 
both filters, indicating isotropic illumination by a central 
source. The maximum per cent polarization is 15 ± 6 and 
10±4 at J- and if -bands respectively (Table|2l • These values 
are lower limits to the intrinsic polarization, since in these 
observations it has not been possible to correct for dilution of 
the polarized flux by the unpolarized light from the central 
star. 

2.2 Hubble Space Telescope observations and 
results 

We have obtained archived HST images for IRAS 
19306+1407^ observed on 2003 September 8 (proposal ID: 
9463). The observations were obtained with the Advanced 
Camera for Surveys (ACS), in conjunction with the High 
Resolution Channel (HRC), using F8I4W- and F606W- 
filters with pivotal wavelengths of 5888 and 8115 A respec- 
tively. The images were reduced using the On-the-Fly Repro- 
cessing of HST Data (otfr) , which produces a cosmic-ray 
cleaned, calibrated, geometrically corrected mosaic image. 
Aperture photometry was performed using GAIA, using the 
Vega zero points^, and obtained magnitudes of 13.81 ± 0.03 
and 12.45 ±0.02 for F606W a.nd F8I4 W respectiYely (Table 

m 

The reduced F606W and F8I4W images are shown in 
Fig. 2 (a) and (b). Fig. 2 (c) shows the F606W image super- 
imposed with contours of J-band polarized flux from Fig. 
1 (b). The object is clearly bipolar in the F606W image, 
and the curved edges of bipolar cavities, extending for 3 to 
4 arcsec from the source, can be seen. The orientation of 
the bipolar axis, at PA 18 deg, is aligned with the J-band 
elongation in total and polarized intensity seen in Fig. 1 (a) 
and (b). The bipolar structure appears to be surrounded by 
a faint, more spherically symmetric halo, seen in both HST 
filters, and this corresponds in extent to the outer contours 
in Fig. 1 (a) and (c). The polarized flux shoulders, at PA 136 
deg, are not perpendicular to the major axis of the nebula 

^ Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the data archive at the Space 
Telescope Institute. STScI is operated by the association of Uni- 
versities for Research in Astronomy, Inc. under the NASA con- 
tract NAS 5-26555. 

^ Ihttp : //www. stsci .edu/hst/acs/analysis/zeropointsl 



Figure 2. High resolution images are available at 
http://star-www.herts.ac.uk/~klowe/ The HST ACS im- 
ages, scaled logarithmically, of IRAS 19306±1407. (a) F606W 
(5888A) scaled between 22 and 13 mag arcsec"^ with the an- 
gle of the major axis indicated by the arrow, (b) F8I4W {8II5A) 
scaled between 22 and 11 mag arcsec"^. (c) F606W ima,ge, scaled 
as above, and J-band polarized flux contours. The lowest contour 
level is 19 mag arcsec"^ and subsequent contours are separated 
by 1 mag arcsec"'^. 



and this is clearly seen in Fig. 2 (c). This non-orthogonality 
in the two axes will be discussed further in Section 4. 

The southern bipolar lobe appears to be the brighter of 
the two in both HST filters, which could indicate that the 
major axis is slightly inclined to the plane of the sky. 



2.3 Sub-millimetre observations and results 

Observations were made on 2005 January 8 using the Sub- 
millimetre Common User Bolometer Array (SCUBA) at the 
15 m James Clerk Maxwell Telescope (JCMT) on Mauna 
Kea, Hawai'i. The SCUBA observations were made simul- 
taneously at 450 and 850 |i.m in photometry mode using 
a jiggle pattern. The 450 and 850 ^m photometry data 
were reduced using the SURF package within the STARLINK 
suite. The sky opacity was corrected using the Caltech Sub- 
millimetre Observatory (CSO) tan relationship*. Flux cali- 
bration was performed using Mars, inclusive of a maximum 
±5 per cent error due to the orientation of Mars' poles rela- 
tive to the Earth and Sun. IRAS 19306±1407 was detected 
at 450- and 850- nm at > la and > Scr respectively inclusive 
of calibration errors. The fluxes obtained (Table for ^450 
and Fsso are 49.9 ± 38.7 mjy and 14.1 ± 3.7 mjy within a 
beam size of 7.5 and 14 arcsec respectively. 



* Using the revised 2000 October 25 relations 
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3 MODELLING THE CSE 
3.1 Model details 

To investigate the dusty CSE around IRAS 1 9306+1407, we 
use modified versions of the iMenardI il989ft axisymmetric 
light scattering (ALS) code to produce Stokes /, Q, U im- 
ages and the axisymmetric radiative t ransfer (DART) code 
iEfstathiou fc Rowan- RobinsonI [l99ol) to model the SED. 
Both codes hav e previously been used to model the CSEs of 
post-AGB stars. iGledhill fc YatesI ll2003h used DART to sim- 
ulate multi- wavelength mid-infrared imaging observations of 
IRAS 17436+5003, in which an axisymmetric shell was re- 
solved. To simulate the axisymmetry, these authors used a 
simple dust density formulation from iKahn fc WestI (Il985l) 
which was found to successfully reproduce all of the axisym- 
metric features, including the offset location of the bright- 
ness peaks seen in the data, which was found to be due 
to the inclinat ion of the system to the plane of the sky. 
iGledhilll (I2OO5I) have used the ALS code to produce generic 
light scattering models of PPN at varying optical depth and 
also find that a Kahn & West density model provides a good 
representation of the observations with a minimum number 
of model parameters. It is important that the dust density 
model uses a minimum number of parameters whilst achiev- 
ing an adaptable axisymmetric geometry, so that there is a 
better chance of each parameter being observationally well 
constrain ed. More complex du st density formulae have been 
used (e.g. iMeixner et 'al]|2002l) . which incorporate the pres- 
ence of AGE and superwind mass loss histories, but require 
more parameters (twice as many in the case of|Moixncr ct al. 
|20o3). These models result in morphologies that are quali- 
tatively similar to our simpler models, but are unlikely to 
be well constrained by our observations. In both the ALS 
and DART models wc therefore use a simpler density profile 
from .Kah n & West (1985) to model an axisymmetric shell, 
whilst recognising its limited ability to reproduce more com- 
plex morphologies: 

pir,9) = po( — Y\l + esm-'e), (1) 

where po is the density at the pole {8 — 0°) at the 
inner radius, rin, and (3 specifies the radial density distribu- 
tion. The azimuthal density distribution is determined by 
parameters e and 7, which specify the equator-to-pole den- 
sity ratio (1 + e) and the degree of equatorial enhancement, 
respectively. An increase to 7 flattens the density distribu- 
tion, creating a more toroidal structure. 

All parameters in Equation 1 are optimized in the 
model, apart from /3, which is fixed at a value of 2 due to 
a limitation of the DART code, corresponding to constant 
mass-loss rate and expansion velocity for the AGE wind. 
The ALS density profile includes an extra parameter, that 
restricts the axisymmetry to within a radius, rsw, modifying 
EquationQto: 

p{r,e) = po (^^"^ when r > r^^. (2) 

A power law size distribution is used with spherical grains 
of radius a, between a minimum and maximum grain size of 
flmin and flmax respectively, and a power-law index, q,: 

n{a) tx a"' for a^in ^ a < Omax. (3) 



The inclination of the symmetry axis to the plane of the 
sky is not known. As mentioned in Section [2.21 the south- 
ern bipolar lobe appears slightly brighter than the north- 
ern one in HST imaging (Fig 2), which could indicate a 
small inclination to the plane of the sky. Although the near- 
infrared images appear consistent with zero inclination (e.g. 
th ey are similar t o edge-on axisymmetric shell models shown 
m 'Gledhill 2005), we consider the inclination angle to be a 
free parameter and allow it to vary in steps of 10 deg. 

The ove rall chemistry o f the system is uncertain. The 
results from iHrivnak et al.l i2000h suggest a C-rich nature 
based on emission features consistent with C-rich PNe. 
iHodge et al.l (|20o3) re-evaluated the mid-infrared spectra 
and classified IRAS 19306+1407 as "UIR features coupled 
with emission from crystalline silicates" suggesting a dual 
chemistry nature. The dust species that have been consid- 
ered in our models are amorphous carbon (amC), silicon 
carbide (SiC) and Ossenkopf cold silicates, and we hav e 
obtained the op tical constants from [ P reibisc h et al.l il993l) , 
iPegQuriJ lll988l) and lOssenkopf et al.l (Il992h respectivelv. 

We ran a total of over 150 ALS and over 300 DART 
models to create a model grid for the free physical parame- 
ters (TableEJ. The minimum and maximum grain sizes were 
investigated from 0.005 to 1 \uxl, with a variable grain size 
spacing typically 0.005 to 0.02 ^.m. The grain size power law 
index was varied between 3.0 to 6.0 at increments of 0.5. 
The radial density fall off exponent is fixed at P — 2 and 
cannot not be varied. The bin widths for the CSE parame- 
ters, common to both models are 1, 2, 10° and 0.1x10"^ for 
the equator-to-pole contrast (e), equatorial density enhance- 
ment (7), inclination angle (6) and the ratio of the inner-to- 
outer radii (nn/rout) respectively. The stellar temperature, 
T, , was investigated using a series of Kurucz models ^ with 
solar metalicities and temperatures separated by 1000 K. 

The ALS code is used to determine the best-fitting enve- 
lope parameters based on the morphology, azimuthal profiles 
in polarized flux and radial profiles of the percentage polar- 
ization and total intensities. The ALS code is additionally 
used to constrain the dust grain size by generating polar- 
ization information. The ALS estimate of the grain size is 
an important input to the DART calculations, which would 
otherwise suffer from a degeneracy between grain size and 
outer CSE radius, both of which strongly influence the long- 
wavelength tail of the SED. The optical depths at 0.55, 1.2 
and 2.2 ^.m are also derived from the ALS model and sub- 
sequently inserted into the DART model. The DART model 
fits to the SED are used to constrain the temperature of the 
central star, inner-to-outer and stellar-to-inner radii ratios. 
The two codes were used to iteratively produce a convergent 
model. 



3.2 Model results 

3.2.1 ALS model 

Before the raw model images can be compared with the 
polarimetric observations, they must be smoothed to mimic 
the effect of the atmosphere and telescope. We find that a 
simple Gaussian filter is unable to reproduce the wings of 

' |http : //kurucz . harvard . edu/ grids ■ html | 
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Figure 3. High resolution images are available at http://star-www.herts.ac.uk/~klowe/ The 1.2- and 2.2-^m smoothed 
model images of IRAS 19306+1407 arc displayed at the top and bottom of the figure respectively. These images are rotated to a PA of 
136°to mimic the observed data. As with the observed images they have been scaled logarithmically. The total intensity (1) is displayed 
in sub-figures (a) and (c) with overlaid polarization vectors (pol) and polarized flux is shown in (b) and (d) . The model images have been 
normalised at the same levels as the observed images: (a) and (c) are scaled between 20 and 13 mag arcsec"^ with lowest outer contour 
levels at 19 and 18 mag arcsec"^, respectively, separated by 1 mag arcsec"^; (b) and (d) are scaled between 20 to 16 mag arcsec"^ and 
19 to 16 mag arcsec"^ respectively with lowest outer contours at 19 (b) and 18 (d) mag arcsec"^, separated by 0.5 mag arcsec"^. 



the PSF effectively, which is essential since the PSF wings 
have a critical effect on the percentage polarization in the 
envelope where the intensity is low, at r > rin. To obtain a 
more realistic fit we use a Moffat filter profile: 

-/3„ot 



M(r) (X 



1 + 



(4) 



where r is radius from the sou rce and ftmof and /3mof are fit- 
ting parameters The Moffat parameters were 
calculated by fitting to the PSF of a bright field star (Ta- 
ble |1J and their uncertainties were estimated by examining 
the fit to the remaining field stars. The filter was then ap- 



plied to the raw (/, Q and U) model images, which were then 
combined to obtain polarized flux and per cent polarization 
values. 

The resulting best-fitting smoothed model is shown in 
Fig. Inland the parameters used are displayed in Table|21 The 
model reproduces the centrosymmetric polarization pattern 
and the observed degrees of polarization in the J- and K- 
bands. The polarized fiux images show the shoulders seen 
in the observations, due to the enhanced scattering at the 
inner edges of the axisymmetric shell, where the dust den- 
sity is greatest. In Fig. 1, the observed polarized flux images 
show a peak of emission at the location of the star. Any mis- 
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Figure 4. High resolution images are available at http: / /star-www. harts. ac.uk/ ~k;lowe/| Left: Azimuthally averaged radial 
profiles of the normalised total intensity. Centre: Azimuthally averaged radial profiles of the per cent polarization. Right: Radially 
averaged azimuthal profiles of the normalised polarized intensity. In all cases, the J- and i^-band data are displayed as squares and 
triangles respectively, with 3cr error bars, and the 1.2- and 2.2-^m smoothed model data are displayed as solid and dashed curves 
respectively. 



alignment of the bright, centrally-peaked images during the 
data reduction stages will lead to a residual polarization at 
this location. Since the polarized flux peak is narrower than 
the seeing disc size, we cannot treat it as significant. We do 
not see polarized emission from the location of the star in the 
model images, since forward-scattered hght (i.e. scattering 
angles close to zero) is strongly depolarized. Higher spatial 
resolution observations will be required to investigate the 
polarization within 0.2 arcsec of the star. If there is signifi- 
cant polarized emission from this region, then an additional 
dust component would be required in the model. 

The fit was assessed by comparing the full grid of ALS 
models to the polarimetric observations. In particular, the 
radial and azimuthal profiles of the smoothed model images 
and the observations were compared and the profiles for the 
best-fit model are shown in Fig. |1| The total intensity image 



radial profile fit (Fig. |1] left) provides a check on the level 
of smoothing, and shows an excellent fit to the observed 
intensity profile at both wavelengths. The fit to the radial 
distribution of per cent polarization (Fig. 2] centre) allows 
us to constrain the dust grain parameters and optical depth. 
The maximum degree of polarization produced by the model 
is very sensitive to the grain size distribution so we consider 
that the grain size is well constrained. The radial distribu- 
tion of per cent polarization depends strongly on the optical 
depth (and hence the dust density), since this determines 
the surface brightness of the CSE relative to the unpolar- 
ized hght from the smoothed PSF. We determine an optical 
depth of 0.68 and 0.11 at J and K respectively, so that the 
CSE is optically thin in the near-infrared. The axisymmetry 
parameters, e and 7 are determined by comparing azimuthal 
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Table 3. The CSE and dust grain parameters for the best-fitting 
ALS and DART models for IRAS 19306+1407. 



Parameter 



Value 



Description 



Dust grain parameters 

Ossenkopf 

Cold Silicates^ 1.0 ± 0.01 
Omin (|Jm) 0.10 ± 0.01 

(nm) 0.40 ± 0.01 

q 3.5 ± 0.5 



Number fraction 
Minimum grain radius 
Maximum grain radius 
Grain size power law index 



Common envelope model parameters 



/3t 
e 

7 

e (deg) 

nn/r-out (10-2) 

ALS model parameters 



2 

6 ± 1 
5 ± 2 
± 10 

7 ± 1 



rU (xlO-i) 



6.78 ± 0.05 
1.13 ± 0.01 
2.0 ± 0.5 



DART model parameters 

T^, (lO^ K) 21 ± 1 

r^Ain (10-5) 1.4 ± 0.2 

A^SE (mag) 2.0 ±0.1 



Radial density fall off 
Equator-to-pole density contrast 
Equatorial density enhancement 
Inclination angle (from equator) 
Inner-to-outer radius ratio 

Optical depth at 1.2 \aa 

Optical depth at 2.2 \aa 

Super-wind 

to inner radius ratio 

Effective Stellar Temperature 
Stellar-to-inner radius ratio 
Equatorial optical extinction 



^Ossenkopf et al.lll992^ . ^This variable is fixed in our model code 
and cannot be varied. tThe optical depth is an output of the ALS 
model. 



Table 4. The Moffat filter profile parameters, Qmof a^nd Anofi for 
a bright field star at J & iC. 



Band 


Ctmof 


/3mof 


J 


3.95 ± 0.06 


2.4 ±0.2 


K 


3.03 ± 0.02 


2.2 ±0.3 



polarized flux profiles to the data (Fig.0]right). The best fit 
gives an equator-to-pole density contrast of 7. 



Table 5. Photometric values f or IRAS 1930 6 +1407 col- 
lated from the l iterature: (1) |H^vna^e^alJ i2000l) : (2) 
Monet et a'i~li200ah: (3) M SX Bands l E gan et al.l2003h . and (4) 
Joint IRAS Science working gr oup 1 1988 



Band 


Central 


Flux density 


Reference 




wavelength 








(vmi) 


(Jy) 




V 


0.55 


7.40 X 10"^ 


(1) 


R 


0.44 


2.21 X 10-2 


(2) 


MSX A 


8.28 


1.16 


(3) 


IRAS 12\Mn 


12.0 


3.58 


(4) 


MSX C 


12.13 


3.65 


(3) 


MSX D 


14.65 


9.12 


(3) 


MSX E 


21.34 


46.27 


(3) 


IRAS 25|jm 


25.0 


58.65 


(4) 


IRASm^aa 


60.0 


31.83 


(4) 


IRAS 100|jm 


100.0 


10.03 


(4) 



Figure 5. High resolution images are available at 
http://star-www.herts.ac.uk/~klowe/ The observed SED 
and best model fits for IRAS 19306+1407. The dash line is the 
model fit and the solid black line is the model fit with interstellar 
reddening a pplied. Refer e nces: (1) thi s paper. (2) [Hrivnak et alj 
1 200( 1. (3) lUeta et al.l ^^i^J£i_^^^^^^^^^^^J^ 
lE garTet alj l20o3r "and (&) Ijoint/^^^^cienc^working groupl 
il988^ . 



3.2.2 DART model 

The SED of IRAS 19306+1407 is plotted in Fig. El using 
published photometry and spectroscopy from a variety of 
sources, including this paper, and covering wavelengths from 
the l/-band through to the sub-millimetre. The photometric 
values are listed in Table |S| The double-peaked nature of 
the SED is immediately evident, consisting of a reddened 
stellar peak around 1.6 \ira and a broad thermal dust peak 
between 30 and 40 (xm due to the CSE. Double-peaked SEDs 
are ty pical of post-AGB stars with optically thin detached 
CSEs Jvan der Veen. Habing fc GeballdlTosa) . 

Our best-fitting model is shown in Fig. 5, both with 
and without correction for interstellar extinction (see be- 
low) . Previous attempts to model the SED using amorphous 
carbon dust and a cooler F/G type st ar, were found not t o 
provide suflicient flux in the dust peak jHrivnak et al.l2000l) . 
We have treated the stellar temperature as a free parame- 
ter and determined a best-fltting value of 21,000 K, typical 
of a Bl type star. This is consistent within erro rs with the 
observationally determined sp ectral type of BO: llVolk et alJ 
120041: iKellv fc Hrivnak l2005l) . where the colon denotes an 
uncertainty in the (Hrivnak, private communication). 

An optical extinction of Ay = 2.0 ± 0.1 mag, through 
the CSE in the equatorial direction, was determined from 
the model fit. We have investigated the effect of inclination 
of the nebula axis and have determined that the SED is 
consistent with a value of 0° ± 10°. The extinction through 
the CSE along our line of sight is, therefore, also Ay ~ 2.0. 

IRAS 19306+1407 lies close to the galactic plane, I = 
50.30° and b = -2.48°, and the SED wiU be aflfected by 
interstellar extinction. The extinction through the Galaxy 
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at this point is estimated to be Ay ~ 5.1 ± 0.2 mag. 
This value was obtained from the IRAS dust reddening 
and extinctio n service^, ba sed on the data and technique 
in .Schlegel. Finkbeiner fc Da vis (199S). 

To correct the emergent model flux for interstellar 
extinction, we apply a reddenin g model developed by 
ICardelh. Clayton fc Mathi3 dlOSSTl which gives the extinc- 
tion, Ax, at every wavelength between 0.1 and 3.3 \im for a 
given Ay and extinction ratio, Ry. The extinction at shorter 
and longer wavelengths has been extrapolated. The DART 
model flux, -Fdart, is then modified to give the flux after 
correction for interstellar extinction. Fx: 

Fx = Fdart X 10 2 = ^ (5) 

Assuming a standard value of Rv=3.1 for the ISM, 
then a fit to the SED shortward of 6 ^.m gives a value 
of 4.2 ±0.1 mag for interstellar extinction (solid curve 
in Fig. 13 . The total extinction to the star is, therefore, 
6.2 ± 0.2 mag. This is consistent with the observed J-K 
colours. Assuming an extinction ratio (Rv) of 3.1, and an 
intrinsic colour excess of E{J — K)o — —0.09 for a BlI star, 
gives Ay = 6.4 ± 0.7 mag. The model parameters used in 
DART are presented in Table |3 

3.2.3 Distance estimate and derived parameters 

The interstellar extinction can be used to estimate the dis- 
tance of the post-AGB star. Using iJoshi (,2005.) . based on ex- 
tinction towards open clusters, gives an estimated extinction 
of 1.58±0.04 mag kpc~^. A visual extinction of 4.2 ± 0.1 mag 
suggests a distance of 2.7 ±0.1 kpc, which we now adopt as 
our assumed distance from this point onwards. 

Using this distance estimate gives values for rin and 
rout of 1.9 ± 0.1x10^* and 2.7 ± 0.1x10^^ m respectively. 
Multiplying rin by r^^/rin gives a stellar radius, of 3.8 ± 
0.6 Rq. 

The stellar luminosity, L*, is obtained by calculating 
the integrated flux under the model SED, giving values of 
1800 ± 140 and 4500 ± 340 Lq, with and without interstel- 
lar reddening applied respectively, for the assumed distance. 
Post-AGB stellar evolution models sug gest a lower limit o f 
2500 Lq for the central star of a PN JSchonberneilll983t) . 
which means that IRAS 19306±1407 must be at least 2.0 kpc 
away to satisfy this criterion. 

To calculate the time scales of mass loss, rin and rout 
are divided by the AGB wind speed. Only the H2 and Ha 
kinematic information are available for IRAS 19306+1407. 
These speeds arise from the shocks and fast winds in the 
post-AGB phase, and are not a true reflection of the AGB 
envelope expansion speed, ther efore we h ave assumed a typ- 
ical speed of 15 km s~^ from .Neri et al.l 11998') . The age of 
the CSE is then 5700 ±160 yrs, became detached 400 ±10 yrs 
and the mass loss lasted 5300 ± 160 yrs. 

The number density of dust grains, A*'o, at rin is cal- 
culated from the optical depth, the extinction cross sec- 
tion of the dust and the CSE thickness. The optical depth 
at 1.2 ^m is 0.678 ± 0.005, giving a value of No = 6.1 ± 
3.0x10"^ m~^. Using A'o and integrating the dust density 
distribution gives the total dust mass (Md), and assuming 

^ |http : //irsa. Ipac . caltech. edu/appllcations/DUST| 



a dust grain bulk density of 3x10"^ kg m ^, gives a value of 
8.9 ±5.0x10"'' Mq. 

The gas-to-dust ratio for this object is unknown and we 
have adopted a value of 200 from Heras fc Hon v ( 2005). The 
total mass of the CSE is then 1.8 ± 1.0x10"^ Mq with an 
average mass- loss rate (M) of 3.4±2.1 x 10"'^ Mq yr'^ The 
derived parameters given in this section are summarized in 
Table il 



4 DISCUSSION 
4.1 CSE geometry 

The polarimetric observations, shown in Fig. 1, have been 
interpreted in terms of an axisymmetric shell with an equa- 
torial density enhancement, which is optically thin in the 
near-infrared. The shell model successfully reproduces the 
observed SED from the V band to the sub-millimetre. As a 
further check on the validity of the model, the ALS code was 
run at the central wavelength of the F606 W fUter to simu- 
late the i/ST observations shown in Fig. 2 (a). The results 
are shown in Fig. 6 and we find that the bipolar structure is 
reproduced, inclusive of the flattened contours in the centre 
of the HST i mage. A single axisym metric shell model, based 
on the simple lKahn fc WestI (Il985ll density distribution, can 
account for the morphology of this object over a wide range 
of wavelengths. The transition from bipolar nebula in the 
optical to limb-brightened shell in the near-IR is due to the 
variation in optical depth through the envelope with wave- 
length. At the wavelength of the i/ST observations, the CSE 
is optically thick along the equatorial direction and so light is 
preferentially funnelled along the polar axes before scatter- 
ing into our line of sight, creating the bipolar lobes. The fact 
that the general appearance and extent of the lobes is repro- 
duced by the model indicates that the density structure of 
the shell, in particular the equator-to-pole density contrast 
of 7, is reasonable. At near-infrared wavelengths, where the 
shell is optically thin along the equator, light is mainly scat- 
tered at the inner boundary in the equatorial plane, where 
the dust density is greatest, creating the shoulders seen in 
polarized flux in our observations. 

Since our model calculations are limited to axisymmet- 
ric geometries, one aspect of the observations that we have 
not been able to account for is the non- orthogonality of 
the polarized flux shoulders, at PA 136 deg, and the ma- 
jor axis of the nebula, at PA 18 deg, illustrated in Fig. 2. 
A similar 'twist' has been detected in the mid-infrared im- 
ages of IRAS 17456±5 003, which has a curving polar axis 
jGledhiU fc Yatesll2003ll . and which was also modelled with 
an axisymmetric dust shell. A further similarity between the 
two objects is the unequal brightness of the polarized flux 
shoulders (see Glcdhill et al. 2001). In the context of our 
model, these are due to scattering at the inner edge of the 
axisymmetric shell, so that the scattering optical depth is 
greater on one side of the shell than the other. Assuming that 
the dust properties are the same throughout the shell, then 
this suggests that there is a greater concentration of dust in 
the brighter shoulder. Further evidence for asymmetric dust 
distributions around post-AG B stars is see n in mid-infrared 
images of IR AS 07134±1005 jDaval et al.iri 998) and IRAS 
21282±5050 jMeixner et ai]ll993ll . ICledhiU fc Yated tOQ^ 
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discuss possible causes for these asymmetries and conclude 
that they may arise due to interaction of the mass-losing 
star with a binary companion, although exactly how this 
hap pens is not c lear. 

IVolk et al.1 (l2004) imaged IRAS 19306+1407 using a 
narrowband H2 filter (2.12 fi.m) and a narrowband A" contin- 
uum filter (2.26 \im), to investigate the molecular hydrogen 
emission. Their continuum subtracted H2 image (their Fig. 2 
left) shows a broken ring with limb-brightened edges, which 
appears cospatial with the central dust structure, at PA 136, 
seen in our polarized flux images. The ring can also be seen 
in their 2.26 \im continuum image, so that they have re- 
solved the dust structure that we see in polarized flux. The 
similarity between the polarized flux and H2 images suggests 
that the scatter ed light and mole cular emission originate in 
the same region. fVolk et al]ll2004tl also detect faint extended 
H2 emission lobes, extending from the ring, corresponding 
to the extended bipolar structure seen in the HST images 
(Fig. oriented PA 18°. It appears that the same axis 
twist seen i n the scattered lig ht images may be present in 
H2 emission. I Volk et al.l^l2004^ suggest that the If 2 ring seen 
in their images coUimates the H2-emitting bipolar lobes. 



4.2 Estimation of the dust mass from our sub-mm 
observations 

The mass of dust in the CSE, Md, can be estimated from 
the IRAS 100 ^m flux, Fioo, and t he SCUBA 850 u.m flux 
We ha ve used the method stated in lGledhill. Bains fc YateJ 
l|2002l) to calculate an estimate of the dust mass from 
our observations. The dust temperature is estimated to be 
146 ± 21 K, using Wien's displacement law, with the peak 
dust emission at 35 ± 5 \ira. The 850 (xm flux value given in 
Tableland _Fioo = 10.03±1.30 Jy, gives an emissivity index 
of 1.3 ± 0.1. The assumed density for a silicate dust grain 
is 3x10^ kg m"^. The total dust mass in the CSE, using 
the assumed distance, is then 4.3 ± 0.7xl0~* Mq, which is 
a factor of ~ 2 less than the value obtained from our ra- 
diative transfer model. The difference may arise from the 
simple assumptions inherent in the sub-millimetre estimate, 
particularly that of an isothermal CSE. The bulk of dust in 
the envelope will be cooler than 146 K (the maximum and 
minimum dust temperatures in the DART model are 130 
and 40 K respectively), and will radiate on the long wave- 
length tail of the SED. An isothermal temperature of 100 K 
would result in a dust mass of 7.2 ± 1.7x10"" Mq. Given 
these approximations, we consider that the two results are 
comparable but that the more rigorous model calculations 
from DART and ALS provide a realistic value for the dust 
mass in the CSE. 



4.3 CSE chemistry 

We have modelled IRAS 19306-1-1407 using a silicate dust 
model, with grain sizes between 0.1 and 0.4 \im, which re- 
produces the shell-like morphology in the near-infrared, the 
observed degrees of polarization and the SED. However, we 
also find that a purely C-rich chemistry (amorphous car- 
bon) using larger grains, typically >0 .6 \im, can reproduce 
the observed polarization jLowe fc Gle dhill 200^ and fit the 
overall shape of the SED, although this produces a poor fit 



Figure 6. High resolution images are available at 
|http:/ /star-www. herts.ac.uk/~klowe/ A comparison of the 
F606W HST image and the raw model image from ALS at the 
central wavelength of the filter. The model image has been 

rotated parallel to the long axis (PA = 18°) to match the HST 
image. The contours are spaced at an interval of 1 mag arcsec^^ 
from the peak value. 



at <1 ^m after interstellar reddening is applied. Amorphous 
carbon also does not reproduce the shape of the SED be- 
tween 10 and 20 \im. We have investigated the possibility 
that silicon carbide could fit the 10-20 \ira region, but find 
that it provides too much flux at 11-12 \im and was in gen- 
eral a poor fit to the SED. These regions are modelled more 
effectively using Ossenkopf cold silicates. 

As mentioned in Section 3.1, the simultaneous pres- 
ence of e mission fro m PAHs and c rystalline silicates 
jHrivnak et al. 2000: Hodge et alj|2004l) suggests that the 
CSE has a mixed chemistry (both O- and C-rich). Our sim- 
ple investigations of mixes of carbon and silicate dust in the 
CSE, show that amorphous carbon significantly dominates 
the SED at less than 1 per cent abundance. This suggests 
that if the 10-20 (xm fits require silicate grains, then they 
must be the dominant dust component. However our mod- 
els do not allow us to segregate the O- and C-rich material to 
have, for example, a region of silicate grains close to the star 
with a largely C-rich outflow at larger radii. Such a config- 
uration has been p roposed to explain o bservations of mixed 
chemistry objects jMolster et al.ll2002fl in which the crys- 
talline emission comes from cool silicates trapped in stable 
circumstellar or circumbinary discs. Matsuura ct al. (200^) 
have shown that in the mixed chemistry post-AGB object 
IRAS 16279-4757 the carbon-rich dust, traced by PAH emis- 
sion, is located in a low-density outflow, while the continuum 
emission is concentrated toward the centre. Although our 
single component model, based on silicate grains, is reason- 
ably successful in reproducing the observations, it is almost 
certain that the chemistry of IRAS 19307-1-1407 involves 
both O- and C-rich material, perhaps spatially segregated 
and with more than one size distribution. 
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Table 6. The derived model parameters at the assumed distance 
of 2.7 kpc obtained from ALS^ and DART* models. 



Parameter 


Value 


Units 


Description 




3.8 ± 0.6 


R0 


Stellar Radius? 


nn 


1.9 ± 0.1 


(lOi-*) m 


Inner Radius^t 




3.8 ± 1.0 


(10") m 


Super-wind Radius^ 


fout 


2.7 ± 0.1 


(10^5) m 


Outer Radius 


LI 


4500 ± 340 


L© 


Stellar Luminosity-t 


No 


6.1 ± 3.0 


(10-3) m-3 


Number density at r^^^ 


Md 


8.9 ± 5.0 


(10-*) Mq 


Total mass of Dust^ 


Av 


4.2 ± 0.1 


mag 


Interstellar extinction? 


^max 


130 ± 30 


K 


Temperature at r?^ 


T ■ 
^ min 


40 ± 20 


K 


Temperature at 



*The apparent luminosity of the star, with applied interstellar 
reddening, is 1800 ± 140 Lq- 



5 CONCLUSION 

We present near-infrared polarimetric images of the dusty 
CSE of IRAS 19306+1407, in conjunction with new sub- 
millimetre photometry and archived HST images. The po- 
larization vectors show a centrosymmetric structure with a 
maximum polarization of 15 ± 6 and 10 ± 4 per cent for J- 
and K-hand respectively. The polarized fiux shows a very 
faint elongated distribution at PA 18° with two bright scat- 
tering shoulders at PA 136° . The object is clearly bipolar in 
archived HST images, with the bipolar axis also at PA 18° 

We model the polarimetric data using an axisymmetric 
light scattering code and a dust model based on sub-micron 
sized silicate grains. The observed polarization features are 
well described by a simple axisymmetric shell geometry, with 
an equator-to-pole density contrast of 7. The same shell 
model is used to fit the SED of IRAS 19306-^1407 from op- 
tical to sub-millimetre wavelengths using an axisymmetric 
radiation transport code, to constrain the stellar tempera- 
ture and radius, the optical depth of the CSE and the mass 
of dust in the CSE. We find that a B-type stellar spectrum, 
with T, =21, 000 K, best describes the SED, confirming pre- 
vious suggestions that the object is an early PN. 

The models give a value for the CSE and interstellar 
extinction of 2.0 ± 0.1 mag and 4.2 ± 0.1 mag respectively. 
We estimate a distance, from the interstellar extinction, of 
2.7 ±0.1 kpc and use this value to derive parameters from 
our models. 

The polarimetric imaging shows deviations from ax- 
isymmetry that are beyond the scope of our model calcu- 
lations. There appears to be a greater concentration of dust 
on one side of the star than the other, plus the axisymmetric 
shell is not aligned with the larger-scale bipolar axis, clearly 
seen in archive HST images. Similar features are seen in 
other post-AGB objects and may result from interaction of 
the mass-losing star with a binary companion. Further evi- 
dence for a binary nature is provided by the probable mixed 
chemistry nature of this object. 
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